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^ . Abstract. Gabor frames for signals over finite Abelian groups, generated by an 

arbitrary lattice within the finite time-frequency plane, are the central topic of 
this paper. Our generic approach covers both multi-dimensional signals as well 
as non-separable lattices, and in fact the multi-window case as well. Our generic 
approach includes most of the fundamental facts about Gabor expansions of finite 
signals for the case of product lattices, as they have been given by Qiu, Wexler- 
Raz or Tolimieri-Orr, Bastiaans and Van-Leest and others. In our presentation 
the spreading representation of linear operators between finite-dimensional Hilbert 
, space as well as a symplectic version of Poisson's summation formula over the fi- 

nite time-frequency plane are essential ingredients. They bring us to the so-called 
Fundamental Identity of Gabor Analysis. In addition, we highlight projective rep- 
resentations of the time-frequency plane and its subgroups and explain the natural 
connection to twisted group algebras. In the finite-dimensional setting discussed in 
this paper these twisted group algebras are just matrix algebras and their struc- 
, ture provides the algebraic framework for the study of the deeper properties of 

' finite-dimensional Gabor frames, independent of the structure theory theorem for 

OO . finite Abelian groups. 
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1. Introduction 



In the last two decades a new branch of time-frequency analysis, called Gabor 
analysis, has found many applications in pure and applied mathematics. The con- 
nection between time-variant systems in communication theory and Gabor analysis 
has turned out to be of great importance for the application of Gabor analysis to 
^ ■ real-world problems such as the transmission of signals between cellular phones. In 

modern digital communication there is an ongoing trend towards FFT-based multi- 
carrier modulation: popular wireline systems for the internet access such as ADSL 
are based on discrete multitone-modulation (DMT) and important wireless systems 
such as WLAN, UMTS, WIMAX make use of OFDM (orthogonal frequency division 
multiplex)-type modulation. In mathematical terminology both DMT and OFDM 
are essentially Weyl-Heisenberg group structured Riesz bases. Thanks to the Wexler- 
Raz relation all the mathematical machinery developed for Gabor frame theory can 
be exploited for the advanced design of so-called pulse-shaped multicarrier systems, 
[¥!] . All these applications lead in a natural manner to the discussion of Gabor 
frames for finite-dimensional Hilbert spaces. Due to the large potential of the mate- 
rial to the above mentioned applications of Gabor frames in signal analysis we want 
to address this note not only to mathematicians but also to engineers who are inter- 
ested in the deeper mathematical background which is likely to provide a sound basis 
for further applications. For the same reason we have tried to present the results in 
a self-contained, although very general form. Finally we want to mention that many 
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researchers have contributed to the construction of frames for finite- dimensional 
Hilbert spaces, e.g. [3JIHI12]. The case of Gabor frames for finite-dimensional Hilbert 
spaces was especially treated by [H El UHl H6l [151 HSl [271 12H1 E21 EH EH HQ] , e.g. the 
Zibulski-Zeevi representation of a Gabor frame operator in the finite-dimensional set- 
ting was discussed in [TO] and in her Ph.D. thesis [31] Ewa Matusiak has presented 
an approach based on representation theory of the Heisenberg group for elementary 
locally compact abelian groups. Recently, Kaiblinger has used finite-dimensional Ga- 
bor frames to approximate dual Gabor windows of a continuous signal, see |25j. 

Let TC be a finite-dimensional Hilbert space of dimension N with inner product 
(.,.)n- A finite family Q = {g;} igI , in TC, with an index set I of cardinality M is 
called a frame if it is a spanning family for TC. Obviously M > N. There are three 
operators naturally associated with Q: the analysis operator Cg, given by 

C g :f ^((f, gi > w ) igI eC M , 

mapping TC into C M and the synthesis operator Dg from C M back TC, given by 

Dg : c = ( Ci ) ie ieC M ^^c igi . 

iel 

An elementary computation shows that 

(Cgf,c) n = (f,D g c) C M , 

hence these two operators are adjoint to each other, Cg = Dg and Dg = Cg. The 
most important operator associated to {gi}iei is the frame operator Sg given by 

Sgf = 

iel 

It is easy to show that Q us a frame if and only if Sg is invertible. In fact, the 
null-space of Sg = Cg o Cg coincides with the kernel of Cg, but ker(Cg) = {0} if 
and only if Q is a frame. Using next the fact that TC is finite-dimensional we can 
argue that the continuous function f i— ► Yliei &)«| 2 is non-zero on the compact 
unit-sphere {f, ||f || = 1} of TC, hence {gijiei is a frame for 7i if (and only if) there 
exist some constants A, B > such that for all non-zero 

Aiifii? f <iic fiij 1I =x;Kf,a>«r<Biifii?,. 

iel 

If A = B in the preceding inequalities, then {gi}i e i is called a tight frame for TC. 
Some authors call tight frames with A = B = 1 normalised tight frames or Parseval 
frames. Any frame Q provides natural expansions of arbitrary elements f G Ti: 

f = s So 1 f = X)< f . s ^a>wa 

iel 

= Sg%f =^(f,gi) W S e 1 g i 

iel 

= S e 1 / 2 SgS, 1 / 2 f = $:(f,S e 1 / 2 g i ) w S e 1 / 2 g i . 

iel 



In the infinite dimensional setting this is the correct definition. 
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These formulas also show how the reconstruction of / from the frame coefficients 
((f) gi)w) ieI is possible. In general the family {gi}; e i is not linear independent which 
implies the non-uniqueness of the frame decompositions of f . 

The paper is organized as follows: In Section 2 we recall some basic facts about 
finite-dimensional matrix algebras. We emphasize the relevance of the existence of 
a trace and of the C*-algebra structure of a finite-dimensional matrix algebra. The 
material in this section makes our exposition self-contained and hopefully accessible 
to engineers and graduate students, even if they are not yet familiar with classi- 
cal Gabor analysis. In Section 3 the finite Heisenberg group is introduced, and we 
present time-frequency shifts as a projective representation of the time-frequency 
plane G x G (analogue to the Schrodinger representation), and its relation to the 
twisted group algebra for G x G. The main result is Theorem 3.8, describing the 
spreading representation for linear operators on C N . In Section 4 we treat (multi- 
window) Gabor systems generated with the help of a subgroup A < G x G and 
characterize the A- invariant operators by their spreading support, respectively by 
their Janssen representations. We also discuss the Wexler-Raz biorthogonality re- 
lations and the duality principle of Ron-Shen for finite-dimensional Gabor frames. 
Finally, Section 5 contains results about the optimality of the canonical dual and 
tight Gabor window with respect to different measures. In particular we point to 
the connection between the tight canonical Gabor window and the Lowdin orthog- 
onalization, [29J. 

2. Basics on matrix algebras 

Let C N and C M be endowed with orthonormal bases {ej : j = 1, ...,N} and {fj : 
j = 1, ...,M}. Then the linear mappings A from C N to <C M can be identified with 
the associated M x N matrix A = (ay) with entries ay = (Aej,f;) for i = 1,...,M 
and j = 1,...,N. We denote with M x N-matrices A^mxn(C) the set of all those 
matrices. The set .Mmxm(C) is an algebra with respect to matrix multiplication, 
with an involution * which for A = (ay) is defined by A* = (a]j) . The map A — > A* is 
an anti-isomorphism, which means that (AA)* = AA* for A G C and (AB)* = B*A*. 
Recall that for general A G A^mxn(C) the matrices AA* and A* A are hermitian, and 
have the same non-zero eigenvalues Ai, ...,A r , counting multiplicities. The square 
roots of the eigenvalues A* A are called the singular values Si > S2 > • ■ ■ > s r > 
where r is the rank of A. A main result on rectangular matrices A in ■Mmxn(C) 
is the singular value decomposition. It tells us that for every A G -Mmxn(C) there 
exist unitary matrices U G .Mmxm(C) and V G A^nxn(C) and a diagonal matrix 
D with non-negative diagonal entries diag(si, s r , 0, 0) such that A = UDV*. 
Despite the possible lack of uniqueness of the factorization, the sequence of singular 
values is uniquely determined, because they are just the eigenvalues of A*A. We 
refer the interested reader to [7] for more information about the SVD. Hence they 
can be used to define a family of norms {||.||sp} on .Mmxn(C), for p G [1, oo], which 
are indeed the finite-dimensional versions of the so-called Schatten-von Neumann 
classes of compact operators on a Hilbert space. They are given by 

i=l 
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The case ||.||s 2 arises in a variety of applications and is called Frobenius, Hilbert- 
Schmidt or Schur norm. In the following we will denote ||A|| S 2 by ||A||r.„ and refer 
to it as the Frobenius norm of A 6 A^mxn(C). The other important Schatten-von 
Neumann norm ||.||sp arises for p = oo. It turns out to be just the operator norm: 

1 1 A 1 1 p = ||A||s°° = lim ||A||sp = si = max{sk : 1 < k < r}. 

p^oo 

The norms || A||p,. and || A|| op are intimately related to its Euclidean structure and the 
C*-algebraic structure of on •Mnxn(C) respectively. These two structures will allow 
us to establish basic facts in a very direct way (also allowing natural counterparts 
in the more elaborate case of continuous variables). 

Recall that an algebra A is a vector space over C with a multiplication compatible 
with the linear structure. If A has a unit element, then A is called unital. If it is 
equipped with an involution A i— > A* satisfying (AB)* = B*A* and A** = A, then 
A is called an involutive algebra. An algebra A is called a normed algebra, if there 
exists a norm on A such that ||AB||_4 < ||A|m|B||„4. We use C*-algebras, which 
are involutive normed algebras (.4, ||.||^) with ||A*A||^ = ||A||^ for all A G A. 

Full matrix algebras A^nxn(C) with the operator norm are C*-algebras, with invo- 
lution given the transition to transposed and conjugate matrices (adjoint operators). 
Conversely, any finite- dimensional C* -algebra A is (iso metrically) isomorphic to a 
full matrix algebra .Mnxn(C), hence finite-dimensionality forces A to be unital. We 
refer the reader to [H] for an extensive discussion of finite-dimensional C*-algebras. 
Next we claim that this are all finite-dimensional C*-algebras. This elementary fact 
will be used in our approach to Ron-Shen duality for Gabor frames in C N below. 

Lemma 2.1. The operator norm is the only operator norm on a full matrix algebra 
turning A^nxn(C) into a C* -algebra. 



Proof. Since A is isomorphic to A^nxn(C) for some N our argument relies on the 
relation between the Schatten-von Neumann classes and the operator norm ||. || op . 
Let A be in A^nxn(C) and D = A*A. We have to show that ||D||_4 = ||D|| op . Let 
di > d2 > • • • d r > be the singular values of A, i.e. the eigenvalues of A*A. The 
equivalence of any two norms on a finite-dimensional vector space yields: 

/ llD 2m || \ / IIDII 2 ™ 

di = inf ( lim — — ) = inf ( lim 4 



R>0 V n^oo R2 m J R> V n->oc R2 m 

In the last equality we have used the C*-algebra property of ||.||,4 and the fact that 
D = D 2 implies ||D 2 |U = hence d x = ||D|U = ||A*A|| op = ||D|| op . □ 



Next, we recall some well-known facts about traces on matrix algebras. For A G 
A4nxn(C) the trace of A = is defined as the sum of its diagonal elements: 

tr(A) = ai ; i H h a NjN . 

The trace tr is a linear functional on A^nxn(C) has the following properties: 



(1) tr(A*) = tr(A). 

(2) tr(AB) = tr(BA) for all A, B G M NxN (C) (tracial property). 
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(1) is obvious, and (2) follows from a simple calculation: 

N 

tr(A*B) = J^ayby = tr(B*A). 

ij=l 

The identification of .Mnxn(C) with C n2 gives that (A, B) i— > tr(AB*) is an inner 
product on .Mnxn(C). In particular tr is non- degenerate, i.e. 

tr(AB) = for all B G M NxN (Q ^ A = 0. 

Therefore (A^nxn( ( C), tr) is a N 2 -dimensional Hilbert space. Consequently, every 
linear functional : A^nxn(C) — > C is of the form 0(A) = tr(AB) for some B G 
A^nxn(C). The tracial property and Schur's lemma, see Proposition 3.6 imply 
uniqueness of the trace. 

Lemma 2.2. Up to multiplicative constants tr is the only linear functional on 
-Mnxn(C) with the tracial property. 

Proof. We have to show that for any linear functional on TWnxn(C) with the tracial 
property 0(AB) = 0(BA), there exists a constant c G C such that 0(A) = ctr(A) 
for all A G -M N xn(C). Let B in _M NxN (C) such that 0(A) = tr(AB*). Then 
tr((BC - CB)A) = 0(AC*) - 0(CA*) = for all A, C G .M NxN (C). Since trace is 
non-degenerate BC = CB for all C G .Mnxn(C), which means that B is in the center 
of A4nxn(C). Consequently B = cIn for some c G C. □ 

Therefore (A, B)p ro := tr(AB*) is the unique inner product on -Mnxn(C) and 

N 

||A||L := tr(AA*) = tr(AA*) = ^ | aiij | 2 

i,j=l 

defines a norm which endows A^nxn(C) with a unique Euclidean structure. 

For S C A the subalgebra generated by the set S is denoted by alg(iS) and is by 
definition the smallest subalgebra of A that contains S, i.e. alg(<S) = spanjsi ■ • ■ s n : 
Si, s n G S}. Recall also that the commutant S' of a set S C TWnxn(C) is 

S' := {A G S : AB = BA for all B G A4 NxN (C)}. 

Observe that S' is always a subalgebra of A^nxn(C) containing all scalar multiples 
of the unit in 7VfNxN(C). Moreover the commutant of .Mnxn(C) is called the center 
nxn(C), see [12] for more information on finite-dimensional algebras. 

Finally, we want to mention some facts about projective group representations. 
Let (G, •) be a finite group and let V be a finite-dimensional Hilbert space. Then 
a projective group representation p : G — > V of G is a family of unitary mappings 
{p(g) : g G G} such that p(g 1 ■ g 2 ) = c G (gi, g2)p(gi)p(g2) for unimodular numbers 
Cc(gi,g2)- The projective group representation p : G — ► V of G is called irreducible 
if {0} and V are the only p-invariant subspaces of V. 
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3. Spreading representation 



In recent years various frame constructions for finite-dimensional Hilbert spaces 
have been proposed, using representations of finite Abelian groups (see e.g. jl]). 
Because these frames inherit symmetries from the groups they tend to have inter- 
esting extra properties. We will study frames associated with representations of the 
Heisenberg group, which is a non-commutative group, a so-called two-step nilpotent 
group, see [HI [39] for a detailed treatment of Heisenberg groups. 

Writing T = {r G C : |t| = 1} for the circle group, the (Weyl)- Heisenberg group 
H(Z N ) = {(r, k, s) : k, s G Z N , r G T} 
is endowed with the following multiplication: 

(1) (n^nXr^Ws) = (riT 2 ■ e 27rik2ri ,k! + k 2 ,n +r 2 ). 



Next we fix some notations concerning Gabor frames for finite-dimensional Hilbert 
space C N . Its elements are considered as discrete time signals of period N, i.e. 
f(k) = f(k + niN) for all m G Z and if k exceeds the period N then k is taken modulo 



N. @ In the sequel we represent a signal f = (f(k)) = (/(0), /(N — 1)) 



as a 

column vector of C N . The Euclidian structure of C N induces an inner product on 
discrete N-periodic signals by (f, g) C N = X]k=o f(k)g(k) for f , g G C N . 

The key players of our investigation are time-frequency shifts of discrete N-periodic 
signals. For an integer k the translation operator T^ is defined by 

T k f = (f(k), f(k + 1), .., f(k - 1)), f = (f(j)) G C N , 

the modulation operator M r is given by 

M r f = (f(0),e 27rir/N f(l),e 2 " i2r/N f(2), e 2 ™ (N " 1)/N f(N - 1)), f = (f(j)) G C N , 

and the time-frequency shift 7r(k, r) of f by 

7r(k,r)f = M r T k f, for f = (f(j)) G C N . 

These operators 7r(k, r) are called time- frequency shift operators (because time shifts 
are applied first) and will be the key-players in our investigation. 

We identify .Mnxn(C) with C n2 and denote the standard basis of C n2 by {<5y : 
i, j = 1, N}. Then the matrix representations of translation Ti and of the modu- 
lation Mi are given by 

/ i o... 0\ ^ 1 o ... \ 

e 2m / N ••• 



Ti 





V i 








1 

• • • J V o ■■• e 27ri(N " 1)/N 

An elementary computation yields T X N = = In and 

MiTi = e 27ri/N T x Mi. 

The preceding equation is a finite-dimensional analogue of the non-commutativity of 
translation and modulation. In his work on the foundations of quantum mechanics 



and M x 



^So in fact C N is identified with € 2 (Z N ) 
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Weyl treated finite-dimensional analogues of the commutation relations which ap- 
peared in the works of Born, Dirac, Heisenberg, Jordan and Schrodinger on quantum 
mechanics, |48j. In these investigations Weyl was heading towards the fundamental 
theorem that all solutions of Born's commutation relation are unitarian equivalent 
to each other. He was only able to prove this result for the finite-dimensional case. 
Later Stone and von Neumann independently turned Weyl's formal arguments into 
a rigorous proof which is the famous Stone-von Neumann theorem. This theorem is 
one of the most important facts in non-commutative harmonic analysis and lies at 
the heart of time-frequency analysis. 

Theorem 3.1 (Weyl). Let U, V be unitary operators on C N such that alg(XJ,V), 
the algebra generated by U and V is all of A^nxn(C), and assume they satisfy the 
commutation relation 

(2) VU = e 2 ^ ik / N UV for gcd(k, N) = 1. 

Then U and V are unitarily equivalent to Ti and Mi, i.e. there exists a unitary 
operator Z such that Z*UZ = Ti and Z*VZ = Mi. 

Weyl had described the result in terms of projective representations of Zn x Zn. 

Proposition 3.2 (Weyl). Let {p(k, s) : k, s G Zn} be an irreducible projective repre- 
sentation of Zn x Zn • Then p is unitarily equivalent to the projective representation 
of x Zn by time-frequency shifts {vr(k, r) = M r Tk : k, r G Zn}. 

It is a fundamental fact of great importance for us that the mapping from (k, r) 
to the time-frequency shift operators {vr(k, r) : k, r G Zn} defines an irreducible 
projective representation of Zn x Zn, i.e. 

(3) 7r(k,r)7r(Ls) = e 27ri(lr - k s) vr(l, s)vr(k, r). 

One advantage of the projective representation instead of unitary representation for 
the Heisenberg group of Zn is that the matrix-coefficients (f , 7r(k, r)g) for f , g G C N 
of {7r(k, r) : 1, r G Zn} have a concrete practical meaning and we do not have to 
take care of additional phase factors. Namely as the Short-time Fourier transform 
of f G C N with respect to a fixed window g G C N 

V B f (k, r) = f 0)g(k - l)e- 27rirl /N = (f , 7r(k, r)g) cN . 
iez N 

An important fact about STFT's is the so-called Moyal Identity. 
Proposition 3.3 (Moyal's Formula). Let fi,gi and f2,g2 in C G '. Then 

(V gl fl, V g2 f 2 ) c | Gx G| = <fl, f2>Cl G l (gl ? g2)cl G l - 

The proof is straightforward and relies on the definition of the STFT and a change 
of summation. Technically speaking Moyal's identity expresses the orthogonality of 
two matrix coefficients of the irreducible projective representation {7r(k, r) : (k, r) G 
G x G} of G x G. An application of Moyal's Identity with f = fi = fi and g = gi = g2 
yields that C = ||g||^| G |- 
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If we regard the elements of Zn X Zn as linearly independent vectors, most con- 
veniently as <5k,o • 5 r ,o f° r k, r e Zn, then their span is the vector space 

A{Z N x Z N ) = | a M)4,o • S r ,o}- 

k,r€Z N 

The vector space ^4(Zn x Zn) has a product which arises from the group product by 
a twist of a uni-modular number. Therefore the arising group algebra is referred to 
as the twisted group algebra of Zn x Zn. More precisely, the multiplication of 5k,(A,o 
and SifiSsfi is the (left) twisted translation 5k,(A,o<5i,(A,o = e27r ^ 1_k ^ N £k+i,(A+s,o- F° r 
the elements of ^4(Zn X Zn) this yields 

( ^ a ( k ' r )^k,0 • 5r,o) ( E S ^ I >° ' 5s '°) 
k,r l,s 

= E E a ( k ' r ) b ^ s)e 2 ™- 1/rm % + i,o • 5 r+s ,o 

k,r l,s 

= E ( E a ^ S W ~ k ' S - r ) e ^ i(1 ~ k>r/ %,o " • 

k,r l,s 

This calculation motivates the definition of a "twisted" product for ^4(Zn x Zn): 

Definition 3.4. The complex vector space A(1i^ X Zn) of functions on the group 
Zn x Zn is called the twisted group algebra when given the twisted convolution \ as 
a product. For a = (a(k, r)) and b = (b(l, s)) in ^4(Zn x Zn) we define the twisted 
convolution of a and b by 

(4) (at]b)(k, r) = E & & s ) b ( ] ~ k ' s ~ r)e 2 " i(1 ~ k) - s . 

l,s 

A projective representation of Zn x Zn induces a representation of the twisted 
group algebra ^4(Z N x Z N ). We define a representation 7r_4 for ^4(Z N x Z N ) by 

Tr^i(a-) := E a (k, r)7r(k, r) for a e *4(Z N x Z N ). 

k,rGZ N 

For non-discrete groups the above representation of the twisted group algebra is 
the so-called integrated representation. The representation tt^ of A(Z^ x Zn) is 
intimately related with the representation of the group Zn x Zn- Therefore we 
collect some of the basic properties of time-frequency shifts which are elementary 
consequences of the commutation relations. 

Lemma 3.5. Let (k, r) and (1, s) be in Zn x Zn- Then 

(i) vr(k,r)* = T_ k M_ r = e 2wik - r / N vr(-k, -r), 

(ii) 7r(k, r)7r(l, s) = e 2 ^- r / N 7r(k + 1, r + s), 
(Hi) 7r(k,r)7r(l,s) = e^^-^^vrO, s)7r(k, r). 

Lemma I3T51 (i) suggests the following natural twisted involution: on A(Z N x Z N ). 

(5) a*(k, r) = e 2 - ik -/ N a (-k,-r). 

Altogether we have the following properties of nx- 
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Proposition 3.6. The mapping a i— > tt a (sl), defined on *4.(Zn x Zn), <?&ven by 

(6) vr^(a) = ^ a(k,r)vr(k,r) 

defines an involutive representation o/.A(Zn x Zn), ie. one has: 

(a) 7r^(a) + n A (b) = tt a (& + b), 

(b) n A (a)n A (b) = Tr^(a^b), 

(c) n A (a*) = ir A (&)\ 

(d) n A (5 fl) = 1 N - 

Conversely, if an involutive representation tt a of A(Zn x Zn) satisfies (a) — (d), 
i/ien t/iere exists a projective representation ir of Zn x Zn satisfying (jHJ), «'.e. ; a// 
suc/i algebra representations arise in this way. 

We still have to show the converse, for which we need the definition of a projective 
representation tc for Zn x Zn by |6|. 

(7) 7r(k,r) := ir A (5 K0 ■ S IfQ ). 

The situation just described is a special case of a more general result, according to 
which there is a one-to-one correspondence between projective representations of a 
finite group G and involutive representations of its twisted group algebra ^4(G x G). 
Another close relation between a group and its group algebra is the fact that tt a is 
an irreducible involutive representation of *4(G x G) if and only if ir is an irreducible 
projective representation of G x G. Schur's lemma formulates this more precisely: 

Proposition 3.7 (Schur's Lemma). Let tt a be an irreducible involutive representa- 
tion o/^4(Zn x Zn) on C N . If a linear mapping A o/C N satisfies 

tt a (a) A = A ^(a) for all a G *4.(Zn x Zn) 

or equivalently 7r(k, r)A = A7r(k, r) for all k, r G Zn- Then A = cIn for some c G C. 

For the proof we refer the reader to the excellent book [IS] by A. Terras, which 
provides further information about representations of finite groups. 

By construction the twisted group algebra ^4(Zn x Zn) coincides with the full 
matrix algebra A^nxn(C). From this point of view Schur's lemma states that the 
center of 7V1nxn(C) consists of all scalar multiples of the identity {cIn : c G C}. 

In our discussion of the spreading function we make use of the existence of a trace 
on .Mnxn(C). The Hilbert space A^nxn(C) with (A, B) Fro = tr(AB*) has a natural 
orthonormal basis which leads to the spreading representation of linear operators on 

<C N 

Lemma 3.8. Let {vr(k, r) : k, r G Z^} be the family of all time-frequency shift 
operators in A^nxn(C). Then {N _1 / 2 7r(k, r) : k, r G Zn} is an orthonormal basis for 

(7W NxN (C), ||.|| F ro). 

Proof. The observation that the cardinality of {vr(k, r) : k, r G Z^} is equal to the 
dimension of -Mnxn(C) and the following calculation gives the desired assertion. Let 
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7r(k, r) and vr(l, s) be two time-frequency shifts for k, r G Zn and (1, s) G Z^. Then 

(7r(k,r),7r(l,s)) = tr(7r(k, r)7r(l, s)*) 

= tr(7r(k - l,r - s)) = N5 k -i, s -r- 

□ 

As a consequence we are able to expand every operator A G .Mnxn(C) with 
respect to all time-frequency shifts from the discrete time-frequency plane Z N x Z N . 

Theorem 3.9 (Spreading representation). For A G .Mnxn(C) we have 

A = (A,7r(k,r)) Fro 7r(k,r) 

k,r€Z N 

= ^A(k,r)7r(k,r), 

k,reZ N 

where i]a = (^A(k, r))k, r gz N called t/ie spreading function of A. Furthermore 

(8) 77 A (k, r) = N- 1 1 " k)e- 27rllr/N /or k, r G Z N . 

iez N 

The expression for the spreading function is a direct consequence of the definitions. 
Traditionally a linear operator A on C N is formulated via a kernel by 

Af(j) = 5> A (i,j)f(i), for f = (f(j)) G C N . 

Then the relation between the spreading representation and the kernel of A is 

j 7A (k,r) = 5]k A (i,i-k) e - 2 ' r<r - i / N 

and the inversion formula is 

kA(i,j) = ^^A(i-j,r)e 2 ^ i / N 

rez N 

The correspondence between A G A^nxn(C) and its spreading coefficients (^(k, r)) 
may be considered as non-commutative analogue of the finite Fourier transform. 
First we have as a non-commutative analogue of Parseval's Theorem which follows 
from the orthogonality of {7r(k, r) : k, r G Z N } and tr(A) = N?7a(0, 0): 

(A, A) Fro = tr(A*A) = N ^ ^(k, r)| 2 for A G A<nxn(C). 

Recall that the Fourier coefficients provide the best least square approximation under 
all trigonometric polynomials. An analogous argument yields the same for general 
orthonormal bases in abstract Hilbert space. For the spreading function this implies: 

Theorem 3.10. Let A be in _Mnxn(C). Then for every subset F C Z N x Z N the 
best approximation to A among all finite linear combinations of time-frequency shifts 
from F - in the Frobenius norm \\-\\fto - ^ given by YIfcz a(k, r)7r(k, r). 
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Proof. The orthogonality of the time-frequency shifts {vr(k, r) : k, r G Zn} implies 
that it is sufficient to prove the following: Let U be a unitary matrix in .Mnxn(C). 
Then A has to be approximated by a scalar multiple of the unitary matrix U. 

|| A - cUHL = tr((A - cU)*(A - cU)) = ||A||| ro - 2 Rectr(AU*) + N|c| 2 

is minimized for c = N _1 tr(AU*). The choice U = vr(k, r) for some k, r G Zn gives 
the desired assertion since rj^k, r) = N" 1 tr(A7r(k, r)*). □ 

The spreading representation was introduced by Kailath in the electrical engineer- 
ing context of time-variant systems [26] completely independent from mathematical 
physics and representation theory. The parallelism with representation theory of 
the Heisenberg group becomes evident through the work of Feichtinger and Kozek 
in [13]. A real world realization of time- variant filters requires a finite-dimensional 
model which can be implemented on a computer. The spreading representation 
allows a symbolic calculus for time-variant filters which in mathematics are called 
pseudo-differential operators. 

Let A be our finite-dimensional model of a time- variant system. Then the spread- 
ing function tja can be considered as a symbol which contains all information about 
the time-frequency concentration of A. Given the spreading representations of 
A = E k ,rez N VaQl, r)vr(k, r) and B = Ek,rez N Vb(K r)vr(k, r), then after an easy direct 
computation of the spreading representation AB reveals that 

(9) AB= Y, (r?A^ B )(k,r)7r(k,r). 

These observations allow us to complement recent work of Wildberger on a sym- 
bolic calculus for finite abelian groups. In [49J the Weyl quantization was considered 
on finite abelian groups and one of the main results states that there is no good 
symbolic calculus for groups of even order. Our results do not rely on the order 
of Zn which is possible because we use a different kind of quantization, namely 
the Kohn-Nirenberg quantization. In [2T] the Kohn-Nirenberg quantization is dis- 
cussed in full generality for locally compact abelian groups, and its relevance for 
time- frequency analysis is shown in [^T] . The presentation given here makes use of 
the finite-dimensionality of the space of time-variant filters, i.e. linear operators on 
C N . 

The singular value decomposition of A of rank r in -Mnxn(C) may be considered 
as the decomposition of A as the sum of r rank-one operators g ; ® hi for g ; , h ; G C N 
for i = 1, r where g ® h denotes the rank-one operator P g ^ 

(10) P g)h f := (g®h)f = (f , h)g for h G C N . 

More concretely, if A has the singular value decomposition UDV* with singular 
values di > d 2 > • • • > d r > 0, then g ; and h; are the i-th column u ; of U and i-th 
row Vi of V, i.e. 

A = diiii <8> vT + • • • + d r u r (8) v^. 
It thus turns out that the spreading representation of A is given by 

r 

A = J2J2 d i (P Ui , Vi ,7T(l,s)) Fro 7T(l,s), 

i=i i,seZ N 
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i.e. the problem is reduced to the spreading representation of a rank-one operator 
Pg,h for g and h as mentioned above. By definition we get 

(11) Pg,h = N- 1 (h,7r(k,r)g) Fro vr(k,r). 

k,rGZ N 

The symbolic calculus for linear operators allows us to transfer properties of the 
operators with relations for their spreading functions. As an example we look at 
the product of two rank-one operators and their spreading representation. Let 
gi, g2, hi, h 2 be in C N . Then an elementary computation gives that 

( 12 ) P gl,hiP g2 ,h 2 = (g2,hi)P gli h 2 . 

The spreading representation of a rank-one operator and the symbolic calculus for 
linear operators yields the so-called reproducing property 

(g 1 ,7r(&,r)hi)^(g 2 ,7r(fc,r)h 2 ) = (g 2 , hi)(gi, ir(k, r)h 2 ). 

In our discussion of Gabor frames we will return to the coefficients of the spreading 
representation of a rank-one operator and place it into the setting of representation 
theory for a pair of twisted group algebras. At the moment we want to emphasize 
the role of projection operators and their associated one- dimensional subspaces. 
Our discussion indicates the connection between our presentation of the spreading 
representation and the work of Calderbank in [22J and Grassmannian frames |42j . 

3.1. Discrete time-frequency plane G x G. So far we have restricted ourselves 
to cyclic groups Zn- This is no real restriction, since the structure theorem for finite 
abelian groups G allows us to move on to the more general setting of a finite abelian 
group G. 

Lemma 3.11 (Decomposition of Finite Abelian Groups). Let G be a finite Abelian 
group of order N. Then G is isomorphic to a direct product of groups Z p N; where 

N = p^ 1 ■ • • p^ k is the prime number decomposition of N, i.e. G = Z p N x x • ■ ■ x Z p N k . 

A character cq of a finite abelian group is a mapping from G — > T, such that 
c g( x + y, w ) — c g( x > ^) c g(Y) u ) f° r ah x ) y ^ G. The set of all characters of G form 
a group with respect to multiplication which we denote by G. A basic result about 
finite abelian groups asserts that G is naturally isomorphic to G. In addition we 
have an explicit knowledge of cg in terms of the building blocks of G, i.e. 

c G (x,u) = c Zl (li,si) ■ • -c Zk (l k ,s k ) x = (b, • • • ,l k ) e G,u = (si, • • ■ ,s k ) G G, 

where we abbreviated Z p N ; by Zi. 

We refer to G x G as the discrete time-frequency plane. We are therefore faced with 
the discussion of twisted group algebras *4(Gi x G 2 ) for Gi for abelian groups Gi and 
G 2 which turns out to be the tensor product A(Gi x Gi)®^4(G 2 x G 2 ). Consequently 

the twisted group algebra A(G x G) = ACL n x x Z nJ ® • • • ® ACL N k x Z N k ) 

Pi p i Pk Pk 

is isomorphic to a product of full matrix algebras M., r ^ri = M. n 1v n x ® ■■■ ® 

| *-J •X v_i I Pi 

A4 N k N k . In the following we denote an element of G x G by (k, r) to emphasize 

Pk x Pk 

that k, r are elements of G = Z n, x • • • x Z n w . 

Pi p k 
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4. Gabor frames 

In this section we discuss Gabor frames from a matrix algebra point of view. If 
G(g, A) is a Gabor frame for C' G ', then the Gabor system G(g, A ) over the adjoint 
subgroup A° is naturally related to the original Gabor system. We will consider the 
twisted group algebras of A and A° and the investigation of their structure allows a 
unified treatment of Gabor frames. The understanding of the twisted group algebras 
for the subgroups A and A° requires some harmonic analysis over the time-frequency 
plane. We therefore develop this prerequisites in the first part of this section. In the 
second part we introduce multi-window Gabor frames and as a motivation we treat 
the set of all time-frequency shifts as the Gabor system Q(g, G x G) and state the 
Resolution of Identity as a " continuous" analogue for reconstruction formulas for a 
general Gabor system. 

4.1. Harmonic Analysis over G x G. Harmonic analysis over the time- frequency 
plane G x G is the study of the Fourier transform of G x G and its properties 
which differ from the Euclidean Fourier transform of G x G. The difference between 
G x G and G x G arises from the symplectic structure of the time-frequency plane. 
We motivate our investigation of G x G with the commutation relations for time- 
frequency shifts: 

tt(L s) [(k, r)] = cg ((k, r), (1, s))c G ((l, s), (k, r))7r(l, s) = ^KM,M)/\Q\^ s ), 

where f2:GxG^GxG denotes the standard symplectic form 

Q((k,r),(l,s)) := 1-r-k-s. 

Therefore the commutation relations for time-frequency shifts endow the discrete 
time-frequency plane G x G with an intrinsic symplectic structure. 

From the algebraic point of view G x G is isomorphic to G x G but from a harmonic 
analysis point of view these two groups are different objects. Namely, the characters 
of G x G are not of the form c Gx g = c g • c g , i.e. the character group of G x G is 
not G x G. Observe that in the time-frequency plane the character group G and G 
are orthogonal to each other, i.e. a rotation by 7r/2 moves G onto G. These facts 
should convince you that the "correct" characters of G x G are Euclidean characters 
rotated by n/2. More precisely, we define a symplectic character c* Gx g| ((k, r), (1, s)) 

for a fixed (1, s) by 

c* GxG| ((k, r), (1, s)) = c G ((k, r), (1, s))c G ((1, s), (k, r)) . 

Lemma 4.1. The character group ofG x G is { c ^ Gx Q| ((k, r), (1, s)) : (k, r) GGxG}, 
i.e. it is isomorphic to G x G. 

This observation motivates the following "symplectic" analogues of translation 
and modulation operators which are actually the translation and modulation oper- 
ators for G x G. Let F be in C' GxG l Then we define the symplectic translation 
operator by 

T^ v) F(k, r) = F((k - k', r - r')) for (k', r') 
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and the symplectic modulation operator by 

M^ 1>s) F((k,r)) = e 2^((M,(i,s))/|G| F ^ k)r ^ 

for (k, r), (1, s) G G x G. As for translation and modulation operators we have the 
non-commutativity of symplectic translations and modulations: 

In other words {T? kr )M/ ls x : (k, r)(l,s) G G x G} is a projective representation of 

G x G x G x G. Therefore we may define the symplectic analogue of the short-time 
Fourier transform and analyze objects on G x G with respect to this symplectic 
short-time Fourier transform. Grochenig uses this object implicitly in his discussion 
of time- frequency localization operators and pseudo-differential operators, [5]. 

The structure of the character group of G x G indicates that the " correct" Fourier 
transform in time-frequency analysis is not the standard Euclidean Fourier transform 
but the symplectic Fourier transform of F G C' GxG ': 

^.(k,r) = |G x Gr 1 / 2 ^c^ Gx g | ((k,r),(l,s))F((l,s)). 

An elementary computation establishes the following property of the symplectic 
Fourier transform. 

Lemma 4.2. The symplectic Fourier transform F s is a self-inverse mapping of order 
two on C' GxG ', i.e. 

_ r -r-2 _ t 

— Si — 1 £\GxG\- 

Since the symplectic characters c ^ Gx q| are obtained from the characters C|qxG| by 
a rotation of ir/2 the Euclidean Fourier transform T and the symplectic Fourier 
transform T s are related by 

^ s = ^oJforJ=(_ I ° G| ^l), 

where J describes the rotation by tt/2 of the time-frequency plane G x G. We want 
to have a Poisson summation formula for the symplectic Fourier transform but what 
is the symplectic analogue of the dual subgroup A 1 - of a subgroup A? After our 
discussion of symplectic characters it should be clear that the correct choice is JA 1 - 
because this set consists of all points (k, 1) such that C|g|(A, (k, r))cjc[((k, r), A) = 1 for 

all A G A. First this implies that the set of all points (k, r) G G x G which satisfies 
this condition is another subgroup A of G x G. In Gabor analysis A° is called 
the adjoint subgroup of A and its relevance was first discovered by Feichtinger and 
Kozek in [14] . Their approach provided an explanation of Janssen's representation 
of the frame operator for separable lattices in M d x R d in [HI [23] and for certain non- 
separable lattices due to Ron and Shen in [37J EE] • For different reasons Rieffel was 
led to consider A , which he calls the orthogonal subgroup of A [36], and recently 
Digernes and Varadarajan came across this object A° and they refer to it as the 
polar of A, [8]. 
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Theorem 4.3 (Poisson summation formula). For any subgroup A < G x G one has: 

(13) ^F(A) = (A)" 1 ^ F(A°) for all F G C |GxG| . 

AeA a°gA° 

An important consequence of the Poisson summation formula for the symplectic 
Fourier transform is the Fundamental Identity of Gabor Analysis(FlGA) . The fol- 
lowing lemma about the symplectic Fourier transform of two STFT's appears at 
different places in the engineering and mathematical literature, see [11] in various 
degrees of generality. Following [43] we call it the Sussman identity. Therefore the 
next result is a finite analogue of Sussman's identity. 

Proposition 4.4. Let fi,f 2 ,gi,g2 be in C' G '. Then 

^[V gl fiV^]((l ) s)) = V glg2 V^((l,s)). 

Proof. 

•^[V gl fiV^]((l,s)) = ^V gl f 1 (k ) r)V g2 f 2 (k,r)e 2 ^« 1 ' s )'^) 

k,r 

= 5>(k, r)f 1; tt(1, s)vr(k, r) gl )(f 2 , vr(k, r ) & ) e ™n((i,s),(k,r)) 

k,r 

= S)f 1; 7T(k, r)7T(l, S) gl ) (f 2 , 7T(k, r)g 2 ) 

k,r 

= (f 2 ,7r(l,s)fi)(g 2 ,vr(l,s)gi). 
In the last line we have applied Moyal's identity 13.31 □ 

The next theorem is just an application of the symplectic Poisson summation 
formula to Sussman's Identity, which is the above mentioned Fundamental Identity 
of Gabor analysis. 

Theorem 4.5 (FIGA). Let A be a subgroup of G x G and fi, f 2 , gi, g 2 G C |G| . Then 

X> Bl fi(A)V^fT(A) = | A |-i V glg2 (A°)V^(A°). 

AeA A°eA° 

In the continuous case Janssen was the first to call this identity the FIGA because 
important results such as the Wexler-Raz biorthogonality conditions and the Ron- 
Shen duality principle are easily derived from it, [23J. 

The relevance of the symplectic structure of the time-frequency plane G x G 
was first pointed out by Feichtinger and Kozek in the context of quantization of 
operators and Gabor frames for elementary locally compact abelian groups j!4j . 
Their approach invokes the full power of abstract harmonic analysis and is therefore 
very technical and not accessible to the majority of workers in Gabor analysis. On 
the other hand these groups are applying Gabor analysis to real world problems and 
one reason for this article is to transfer the main results of Feichtinger and Kozek 
into a form which requires only a modest background in harmonic analysis. 
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4.2. Multi-window Gabor frames. In this section we explore mult i- window Ga- 
bor frames with the help of the spreading representation. This allows for a uni- 
fied discussion of non-separable Gabor systems which generalize the known results 
for separable Gabor systems due to Wexler-Raz, Tolimieri, Qiu and Strohmer, 
[171 1351 0:0] • First of all we treat just Gabor frames and obtain their main prop- 
erties and at the end of this section we indicate how these results are extended to 
multi-window Gabor frames, i.e. a finite sum of Gabor frames. 

Definition 4.6. Let A be a subgroup of G x G and g a Gabor atom in C' G '. Then 
£?(g, A) = {7r(A)g : A G A} is called a Gabor system. If the frame operator 

S giA f = ^(f,vr(A)g)7r(A)g 

AeA 

is invertible, then £?(g, A) is called a Gabor frame. 

The signal g := S gA g is the canonical dual window and h := S gA g is the 
canonical tight window of the Gabor frame Q(g, A). 

We consider the elements of the Gabor system {7r(A)g : A G A} as columns of 
a |G| x |A|-matrix D gi A = [7r(Ai)g, 7r(A|A|)g] for some ordering of the elements 
of A. Then the Gabor frame operator S g; A may be written as S gi A = D gj A o D* A , 
i.e. S gi A acts on vectors in C' G ' and implements the coefficient operator for Q(g, A). 
Observe that the operator D g A o D g A with respect to the canonical basis {5a,o : 
A G A} is represented by the Gram matrix G gi A of Q(g,A) with entries (G^ iAt = 
{&[f> §a)) ^ g A - Consequently the Gram matrix G g A acts on vectors in C' A L 

As a motivation of a general Gabor system Q(g, A) we first treat the Gabor frame 
£?(g, 0}), i.e. {vr(k, r) : (k, r) G G x G}. Then the coefficient operator and synthesis 
operator are given by 



Cgf = ((f,7r(k,r)g)) andD g c= V c(k, r)7r(k, r)g for c G C |G> ( ;| 

V / (k,r) 



(k,r)GGxG 



Proposition 4.7. For g G C' G ' we have that Q = {vr(k, r)g : (k, r) G G x G} is a 

tight frame for C' G ' with frame constants A, B = ||g||^| G |, i.e. 



fllcic = E Kf,^(k,r)g>| 2 for all f G C' G I. 



|2 

I^IGI,. , 

G x G 



k,rgG 



Corollary 4.8 (Resolution of Identity). Let g G C' G ' with ||g||ciGi — 1- Then for 
every h G C' G ' with (g, h) 7^ one has 

f= |G|.L h ) E (f^(k,r)g)vr(k,r)h. 

(k.r)eGxG 

The corollary follows from our discussion of frames in the introduction with ho : = 
Sg g. Our proof that {vr(k, r)g : (k, r) G G x G} is a tight frame for C' G ' relies on 
the commutation relations for time-frequency shifts and Schur's lemma. The frame 
operator Sg of {7r(k, r)g : (k, r) G G x G} has the following property: 
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Lemma 4.9. For every (1, s) in G x G we have 

Sg = 7r(l,s) oSgO 7r(l, s)*. 

Proof. The commutation relations 7r(l, s)*7r(k, r) = e _27r ^ k_1 ^ s /l G l7r(k — 1, r — s) yields 
7r(l,s)oS o7r(l,s)*f = ^(7r(l,s)*f,7r(k,r)g)7r(k,r)g 

k,reG 

= Yl < f ' < k ~ *> r - 8 )sM k - 1, r - s)g = S e f . 

k,r€G 

□ 

Recall that {7r(k, r) : (k, r) : G x G} is an irreducible representation of G x G 
and that it generates A4\g x g\(C), hence the preceding observation yields that Sg is 
in the commutant of .M|GxG|(C). Hence by Schur's lemma Sg is a multiple of the 
identity operator, 

Sg = C ■ I|g| for some C G C. 
The determination of the constant C follows from an application of Moyal's identity 
This implies the assertion that {7r(k, r) : (k, r) G G x G} is a tight frame for C' G L 

Before we move on to general Gabor frames we note some important properties 
of the spreading representation, which allows us to justify the name "Resolution of 
Identity" . First of all we remark that the projective representation of G x G gives 
rise to a unitary representation on -M|gxG|(C). In other words 

(k,r) h- A[(k,r)] = 7r(k,r)A7r(k,r)*, A G M| G xG|(C) 

is an involutive automorphism of A4\ GxG \(C). 

Proposition 4.10. The mapping (k, r) i— > A[(k, r)] defines a unitary representation 
of G x G on the Hilbert space .M| GxG |(C) with the Frobenius norm, i.e. 

(1) A[(k,r)]oA[(l,s)]=A[k + l,r + s] ; 

(2) (A[(k,r)],B[(k,r)]) Fro = (A,B) Fro . 

We are interested in the relation between the spreading function of an operator 
A G .M|GxG|(C) and of A[(k, r)] for (k, r) G GxG. Since the spreading representation 
of A is an expansion with respect to the basis {vr(k, r) : (k, r) G G x G} the problem 
reduces to the understanding of conjugation by 7r(l, s) for time-frequency shifts. 

Lemma 4.11. For A G M.\gxG\(^) and (k, r) G G x G one has ?7A[(k,r)] = M^t^a- 
Proof. 

7r(k,r) o Ao7r(k,r)* = ^ VaQ, s)7r(k, r)?r(l, s)?r(k, r)* 

l,s 

= ^ ((k, r), (1, s))x G ((l, s), (k, r))r ?A (l, s)vr(l, s) 

l,s 

= E M (1,S)^(1,S). 

1,8 
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□ 

This behaviour of the spreading coefficients under conjugation of the operator will 
be crucial in our discussion of the Janssen representation of Gabor frame operators. 
At the moment we want to point out that the preceding properties of the spreading 
representation gives the following form of the Resolution of Identity ( 14.81) . Namely, 
recall the rank-one operator P gi hf — (f>g)h for f , g G C' G l Then an elementary 
calculation gives that 

[vr(k, r) o P g>h o 7r(k, r)*] f = (f , vr(k, r)g>7r(k, r)h. 

Therefore conjugation of a rank-one operator P g h by a time-frequency shift moves 
it to the point (k, r) in the time-frequency plane C' GxG L We denote conjugation by 
7i"(k, r) of a linear operator A on C' G ' by 

A[(k,r)] := 7r(k,r) o Ao7r(k,r)*. 

Consequently the Resolution of Identity can be expressed in the following way, 

' k,r 

for (f , h) ^ 0. If g = h then the identity operator is a linear combination of 
orthogonal projections onto the one- dimensional spaces generated by the family 
{7r(k, r)g}. Recall that the pure states of -M|gxG|(C) are the rank one operators. 
Therefore the resolution of identity may be understood as shifting a pure state of 
■M|gxg|(C) over the discrete time-frequency plane C' GxG L 

After these preparations we want to explore the structure of Gabor frames Q(g, A) 
for A < G x G, i.e. discrete analogues of the Resolution of Identity: 

f = ^ (f, vr(A)h) 7r(A)g for a suitable h G C |G| . 
aga 

The last equation indicates that operators of the type 

S g ,h,Af = X)<f,7r(A)h)7r(A)g 

AGA 

are closely related with the frame operator of a Gabor system Q(g, A). Due to 
this fact they are called Gabor frame-type operators. The following property of 
Gabor frame- type operators is crucial for an understanding of the structure of Gabor 
systems and is the very reason for all duality principles in Gabor analysis. 

Lemma 4.12. Let A be a subgroup of G x G and g a Gabor atom in C' G '. Then 
the Gabor frame-type operator S gi h,A commutes with 7r(A) for all A G A, i.e. 



(14) 



7r(A) O S gjh ,A O 7r(A)* = S gj h,A- 
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Proof. The commutation relation for time- frequency shifts yields for fi G A.: 

M^W^)*]* = ^(7r( / i)^7r(A)h)7r( / ,)7r(A)g 

AeA 

= ^(f,7r( M )7r(A)h)7r( / ,)7r(A)g 

AeA 

= ^(f,7T(/i + A)h)7r(/i + A)g = S g , hiA f. 

AeA 

□ 

Since < HM is a crucial property of S g> h,A we want to explore the structure of the 
set 13(A) of all linear operators A G M\gxg\ which commute with all time-frequency 
shifts 7r(A) from a given subgroup A of G x G. We call these operators A-invariant. 

Proposition 4.13. For every subgroup A < G x G £/ie se£ 23(A) zs an involutive 
subalgebra o/M| GxG |(C). 

Proof. If A, B are A-invariant operators, then we have that: 

(1) A + B G B(A) and cA G 13(A) for all c G C 

(2) AB G B(A) since vr(A)A7r(A)* = A and 7r(A)Bvr(A)* = B implies that 

AB = 7r(A)A7r(A)*7r(A)B7r(A)* = 7r(A)AB7r(A)*. 

(3) A G B(A) implies 7r(A)A*7r(A)* = A*. 

□ 

The algebra 13(A) of A-invariant operator A is the commutant of the twisted group 
algebra A(A) within A4\q\ is equal to the twisted group algebra ^4(A°). In detail: 

Proposition 4.14. Let A be a subgroup of G x G. Then we have 

(1) B(A) = A(A°); 

(2) The commutant of 13(A) is A(A) and the commutant of A(A°) is A(A); 

(3) The center of 13(A) = A(A n A°); 

(4) 13(A) is commutative if and only if A° C A. 

The last assertions indicate that the structure of A is essential for the properties 
of 13(A). In analogy to symplectic vector spaces we call the subgroup A of G x G 
isotropic if the symplectic form Q vanishes identically on A. The largest subgroup 
of G x G with this property is naturally called maximal isotropic. A moment of 
reflection shows that a A is isotropic if and only if A° C A. Therefore 13(A) is 
commutative if and only if A is isotropic. Now the maximal commutative subalgebra 
of Ai\GxG\(^) is the algebra of diagonal matrices which implies that A is maximal 
isotropic if and only if A is a product lattice A x A 1 - for A < G and A 1 - < G. In other 
words, the Gabor frame operator S g AxA ± for a product lattice A x A 1 - is unitarily 
equivalent to a diagonal operator. The Zak transform is the interwinding operator 
which diagonalizes the Gabor frame operator Sg^xA- 1 -) see [IS]- 
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Note that : A G A} defines a reducible projective representation of A if A is 

a proper subgroup, because then its commutant ^4(A°) is non-trivial. 

Proposition 4.15. Let A be a subgroup o/GxG. Then {n(X) : A G A} defines a 
reducible projective representation of A. 

The previous statement is valid for any subgroup A of G x G, especially for A°. 
The twisted group algebras »4(A) and ^4(A°) are the matrix algebras underlying 
Gabor analysis over finite abelian groups. Therefore we investigate their structure 
in detail. We use the notion of matrix coefficients for a twisted group algebra. Let 
^A(A) be a representation of the twisted group algebra A. Then we define a matrix 
coefficient of A(A) as 

(7U (A) (a)g,h) d A| =5^a(A)(7r(A)g,h) c |A| for a = (a(A)), g, h G C |A| . 
AeA 

There is a close relation between the matrix coefficients of A(A) and ^4(A°). 
Theorem 4.16. Let A be a subgroup o/GxG. Then we have 

n n 

^^d i (g i , 7 r(A)h i )( 7 r(A)g,h) = — E d i<S^( A °)s>M A °) h ' h i>- 

AeA i=l ' ' A°eA° i=l 

Proof. Recall that every A G A(A) may be written as 

n 

k A (a)(Va) = ^^d i (g i ,7r(A)h i )7r(A) 

AeA i=i 

and this yields to 

n 

(7U(A)(??A)g,h) c |A| = ^^d i (g i ,7r(A)h i )(7r(A)g,h) 

AeA i=i 

I I A°eA° i=l 

where we applied in the last equation the FIGA. The right side of the last equation 
may be understood as the matrix coefficient for a certain element of 7T4(a<>). □ 

If we consider the special case A = G x G, then the statement of the last theorem 
specializes to Moyal's Formula, i.e. the Schur orthogonality relations for STFT's. In 
this sense we consider the last theorem as a generalization of Schur's orthogonality 
relation to reducible group representations. 

The preceding observations yield the following representation of a A-invariant 
operator, which includes the representation of Gabor frame operators as given by 
Tolimieri-Orr, Qiu, Strohmer, and Wexler-Raz for product lattices Ai x A 2 . 

Theorem 4.17 (Janssen representation). Let A be a subgroup of G x G. Then for 
a A-invariant operator A we have a prototype matrix P such that 

A = E p [ A ] = E 7r ( A ) ,toPo7r ( A ) 

AeA AeA 
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is the A-periodization of P or 

A= J2 (P,vr(A°)) Fro 7r(A°), 

A°eA° 

i.e. the spreading coefficients of A are the sampled spreading coefficients of the 
prototype operator P to the adjoint lattice A° . 

Proof. An application of the commutation relations for time-frequency shifts in the 
following form 

TT(fi)n(X)n(^)* = c? Gx g,(A, /u)7r(A — /i) for all A, \i G A 

implies that 7r(/i)o7r_4(A)(a)o7r(/i) is another element of A(A) for a suitable translated 
and shifted version of a. In terms of the spreading representation, the imposition 
of A-invariance on a linear operator A implies that the spreading coefficients are 
periodic, i.e. A G £>(A) is equivalent to A [A] = A, i.e. that for all A G A one has: 

f?A[A](» = C|G|(A,/i)cjo[(^, A)t7a(a0 for all /i G A. 
In other words every A-invariant operator has a representation of the form 

A = a(A>(A°) 

A°eA° 

for some vector a = (a(A )) G C |A °L 

The prototype operator P is the sum of projection operators arising from the 
singular value decomposition of A, i.e. P = J2\=i djUi <8> vj". Let A be the Gabor 
frame operator of Q(g,A). Then 

S g , A =|A|- 1 (g ) vr(A°)g)7r(A ). 

A°eA° 

□ 

We interpret the Janssen representation of a A-invariant operator in terms of 
multi-window Gabor frames. At the moment we draw some consequences of the 
Janssen representation which generalize results of Wexler-Raz on Gabor frame oper- 
ators for separable lattices in cyclic groups. In their work on the inversion of Gabor 
frame operators over cyclic groups Wexler and Raz gave the impetus for the duality 
theory of Gabor frames which was developed independently by several groups of 
researchers. Our general results on the structure of A-invariant operators lead natu- 
rally to the study of a finite number of Gabor frames. Let gi, g r G C' G ' be Gabor 
atoms of the Gabor systems <?(gi, A), £/(g r , A), then we speak of a multi-window 
Gabor system. If the associated frame operator 

r r 

j=i j=i AeA 

is invertible, then the system Q = {J T j =1 G{gj, A) is called a multi-window Gabor 
frame. Therefore our discussion of the Janssen representation of a A-invariant op- 
erator A states that A is a multi-window Gabor frame operator with gi = Uj for 
i = 1, r. We refer the reader to [201 ED E] for further information on multi-window 
Gabor frames. 
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The multi-window Gabor frame operator is actually a finite sum of rank-one 
operators 

Se = $^(gi ® gT[A] + • • • + gr <E> &[A]) 

AGA 

and is a A-invariant operator and therefore has a Janssen representation. 

The main problem in Gabor analysis over finite groups is an understanding of 
all dual pairs (g, h) for a given Gabor system Q(g,A). In other words, we look for 
all (g, h) such that S gi h,A = Ig- We attack this problem by an application of the 
Janssen representation of A-invariant operators. 

Corollary 4.18. Let A be a subgroup of G x G and £(g, A) a Gabor system for 
g G C' G '. Then the following are equivalent: 

(1) f = E AeA <f > 7r(A)h>7r(A)g for all f G C' G ' . 

(2) (g, 7r(A°)h> = |A| • 5 A o j0 /or a// A G A°. 

Proof. The equivalence of (1) and (2) follows immediately from the uniqueness of 
the spreading representation, the fact that 77n |G = oa°,0) an d the identity 

S g ,h,Af = |A| _1 J2 (g' 7r (A°)h)7r(A°)f, for all f G C |G| . 
A°eA° 

□ 

The biorthogonality relations of Wexler-Raz are a relation between traces of -4(A) 
and v4.(A°). Namely, the commutation relations for time- frequency shifts and the 
tracial property yields that 

tr4 (A) (A) = 77a(0) = a , and tr^ (A o } (A) = ci] A (0) for some c G C. 

The Poisson summation formula for the symplectic Fourier transform allows us to 
identity c with 1/|A|, i.e. the traces of -4(A) and A(A°) are multiples of each other: 

1 

tiU(A) = j-^-j- tr^(A-) • 

Let A be a Gabor frame-type operator. Then the preceding equation yields the 
Wexler-Raz biorthogonality relations. More generally, we get a "weighted" Wexler- 
Raz biorthogonality relation for multi-window Gabor frame-type operators. 

The Janssen representation of a Gabor frame operator provides an alternative 
route to the inversion of the Gabor frame operator: Namely, a Gabor frame operator 
may be considered as a twisted convolution operator on A°. More generally, 

Lemma 4.19. Let A be a subgroup of G x G and A, B G 13(A)- Then the spreading 
coefficients of AB are given by the twisted convolution of those for A and B; 

WA°)= ^a(A -/x>b(^ )cg(A -/x ,/x°). 

^°eA° 
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Therefore we have a symbolic calculus for A-invariant operators. We have that 
S g ,A xI A = S gj A, i-e. 

Vs sA (\°) = Vs s Ja°Vi a (^°) = |Ar 1 ^(g,7r(A -/i°)g)c G (A -/i , / i ). 

On the level of representation coefficients the discussion of Gabor frames corresponds 
to the study of the mapping 

c = (ca°) i— > Gc 

with G = (G(A°, /i° G A )) = ((g, tt(A° - /x°)g)c G (A° - fi°, . By the definition of 
the Gram matrix and the commutation relation, G is the Gram matrix of {7r(A°)g : 
A° G A°}. Now the question of invertibility of the Gabor frame operator translates 
into the invertibility of the twisted convolution of the Gram matrix of C?(g, A°). 
More precisely, Q(g, A) is a Gabor frame for C' G ' if and only if G is invertible. Let 
c be a vector in C' A °L Then c represents the canonical dual window ho = S gA g, if 
it is the unique solution of the system of equations 

Gc = |A|<? A o )0 . 

Consider the algebra of all A-invariant operators 13(A) as a C*-algebra of matrices, 
i.e. we equip .M|g|x|A| with the operator norm. Then we may express an element 
of 13(A) in terms of {7r(A) : A G A} or {ir(\°) : A° G A°}. Now by the spreading 
representation we obtain a relation between the operator norm of || S s ,a ||op and ||G|| op . 
By the uniqueness of the C*-algebra norm we get the existence of a constant k such 
that ||S gj A||op — k 1 1 G 1 1 op which for ||S gj A||op equal to the identity yields k = |A| _1 . 
As a summary of the previous statements we get the other duality result for Gabor 
frames, the Ron-Shen duality principle. 

Theorem 4.20 (Ron-Shen duality). Let A be a subgroup ofGxG. Then Q(g,A) 
is a frame for C' G ' if and only if G(g, A°) is a (Riesz) basis for C' A °'. Let Aa,B\ 
be the frame bounds ofQ(g,A) and Aa°,-Ba° the (Riesz) basis bounds of Q(g, A ). 
Then 

A\o = \A\A\ and B^o = \A\B^. 

5. Dual windows, canonical tight windows and Lowdin 

orthogonalization 

In this section we want to stress the importance of the canonical dual and canonical 
tight Gabor window, especially the relation between the structure of the canonical 
tight Gabor frame and the Lowdin orthogonalization of the original Gabor Riesz 
basis. 

Let G(g, A) be a Gabor frame for C' G L Then we denote by r g A the set of all dual 
pairs (g, h), i.e. 

Tg,A = {h : C g A C g)A = I|G|}- 
The exploration of this set is one of the core problems in Gabor analysis. Observe 
that the A-invariance has as important consequence, the A-invariance of arbitrary 
powers of the Gabor frame operator: 

(15) tt(A) o S^ a o n(\y = Sl A for all a G K. 
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Therefore 

f = S g , A S g , A f = £<f, 7 r(A)g) 7 r(A)(S gi 1 A g) 
AeA 

and £?(g, A) is called the canonical dual Gabor frame given by g = S~\g. The duality 
theory of Gabor frames implies that 

S ^ = m £<&*-(A )gMA ), 

i.e. S gi A is an element of the linear span of Q(g,A°) = {7r(A°)g : A° G A°}. By the 
Wexler-Raz biorthogonality relation we can identify r g ,A- 

Theorem 5.1 (Wexler-Raz). Let G(g, A) be a Gabor frame for C' G ' . Then h G r g) A 
if and only if 

Proof. Our argument closely follows |20] . Every h G r g) A satisfies the Wexler-Raz 
biorthogonality conditions, therefore we have (h — g, 7r(A°)g) = for all A° G A°, 
i.e. h — g G Q(g,A°) ± . Conversely, letting h G g + Q(g, A )- 1 , then h fulfills the 
Wexler-Raz biorthogonality relation, i.e. S gi h,A = I|g|- D 

The canonical dual Gabor atom g = S g A g is a distinguished element of r g A since 
it is the "optimal" dual window in various ways. 

Theorem 5.2. Let Q(g, A) be a Gabor frame for C' G '. Then for all h G r g>A the 
following (equivalent, characteristic properties) hold true: 

||g|lci G i — l|h||d G i {minimal norm), 

IIQ|,Af ||ci Q i — ll^h a 1 " llci G i (minimal norm coefficients) , 

||g — g|| c ici < ||g — h|| c |G| (closest to Gabor atom), 

l|g/l|g|lciGi - g/||g||ciGi|| c iGi < ||g/||g|| c iGi - h/ll h llciGi|| c iGi (most likely). 

We refer the interested reader to [33] for the proof of the preceding theorem and 
for a detailed discussion of the structure of r g A and that it actually characterizes 
the dual Gabor atom. We only want to focus on the minimal norm property of the 
coefficient mapping for the canonical dual Gabor frame. The following computation 
shows the relevance of the Moore-Penrose inverse for Gabor frames. 

C g)A f = ((f,vr(A)S gi 1 A g)) = ((f,S gi 1 A 7r(A)g)) 

= «s g)A f, Tr(A)g)) = c kA (q A c^)-H = C g + A f , 

where C^ A denotes the Moore-Penrose inverse of the coefficient operator C gi A. Then 
it is well-known that C^ A is the minimal solution of the system C giA f = a for a given 
coefficient vector a and such that f may be expressed in terms of ((f , 7r(A)g)), i.e. 
f = Ba under the constraint that 

||C g)A f - a||J,| | = ||(C g)A B-I| G |)a||^ G | = minimal. 

In other words we have to minimize ||C g)A B— Iiqi ||| ro over all B which is solved by the 
Moore-Penrose inverse B = C g A . In [T7| it has been shown that the Moore-Penrose 
inverse minimizes 

||C g)A B - J|| Sp 
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for all Schatten-von Neumann norms, ||.||s p for all p G [l,oo]. Consequently, the 
coefficients of the canonical dual Gabor frame are minimal for all £ p -norms. 

There exists a canonical way to associate a tight Gabor frame with Q(g,A). By 
the A-invariance of S gi A we have 

Sg,fS g ,AS^ = ^(f,7r(A)S gi 1 A /2 g)7r(A)S gi fg 

AeA 

Observe that G(h°, A) is a tight frame generated by the canonical tight Gabor atom 

— 1/2 

h° = S g A g. More generally, we investigate the set of all tight Gabor frames Q (h, A) 
for a given subgroup A in C' GxG ', i.e. r tig ht(g, A) = {h 6 C' G ' : Sh,A = I|G|}- 

Theorem 5.3. Assume that h e r«^(g, A). Then there exists a unitary operator 
U such that UCh,A = Ch°,A- 

Proof. The statement of the theorem may be expressed as follows. If CJjo A Ch°,A = 
C£ iA Ch,A holds, then there exists a unitary U such that Ch°,A = UCh°,A- We de- 
fine U on the range ran(Ch°,A) by UCh°,Ai" = Ch.Ai" and by zero on the orthogonal 
complement. Then C£o A Ch°,A = C£ A Ch,A implies ker(Ch°,A) Q ker(Ch,A)- Conse- 
quently U is an isometry satisfying UCh,A = Ch°,A- D 

The previous theorem corresponds to the fact that there is some freedom in taking 
square-roots of a positive operator. The canonical tight window h° is the tight 
frame which minimizes ||g — h|| c |Gi among all tight frames. By the Ron-Shen duality 
principle this is equivalent to the minimization over all orthonormal basis Q (h, A) for 
the linear span of G(g, A). The problem for the orthonormal basis was first solved 
by Lowdin in [29J and corresponds to the choice A = G -1 / 2 = (C gjA C gi A) _1 ^ 2 . 
We formulate Lowdin's orthogonalization for a general set of linearly independent 
vectors Q = {gi,...,gN} of C M . Let C be the M x N matrix with {gi,...,g N } as 
columns. 

Theorem 5.4 (Lowdin). Let Q = {g l5 ...,gN} be a collection of linearly independent 
vectors in C M and singular value decomposition C = UDV* of C. Then ||C — X||F ro 
is minimized for unitary X by L = C(C*C) -1 / 2 = UV*. 

Proof. First we show that L is unitary and secondly we show that it is the unique 
minimizer. Since C = UDV* we have that 

l = udv*(vditudv*;t 1/2 = udv*(vd 2 v*;t 1/2 
= udv^vdv*)- 1 = uv*. 

Since U, V are unitary, we have that L is unitary. 

|| C - X|| Fro = ||UDV* - X|| Fro = ||D - U*XV||fvo = min, 
if U*XV = I, i.e. X = UV*. □ 

In [12] it is demonstrated that L is the minimal orthogonalization for any unitarily 
invariant norm, especially for all Schatten-von Neumann classes ||.||sp andp e [0, oo]. 
Many authors call the Lowdin orthogonalization the symmetric orthogonalization 
because it is invariant under a permutation of the linearly independent vectors of 
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G. More precisely, let P be a N x M permutation matrix. Then LP is the best 
orthogonalization of CP. The proof is by contradiction. 

We close with a few words on the case of Gabor frames. Since the Gram matrix 
of a Gabor frame has a certain structure, the absolute value of entries are circulant, 
and the Frobenius norm becomes the norm of the Gabor atom, see Janssen and 
Strohmer for a more detailed discussion [21] . 

6. Conclusion 

The main goal of this paper is both a self-contained description and a survey 
of Gabor frames for finite-dimensional Hilbert spaces including its duality theory, 
such as the Wexler-Raz biorthogonality relations, Janssen's representation of the 
Gabor frame operator and the Ron-Shen's duality principle. The presentation relies 
on elementary properties of twisted group algebras, which are matrix algebras in 
the given finite-dimensional setting. Our focus on twisted group algebras allows 
an elementary and self-contained approach. The interplay of these twisted group 
algebras underlies the well-known results derived by Wexler-Raz, Qiu, Strohmer and 
others on Gabor frames for product lattices. Therefore we are able to generalize their 
results to arbitrary subgroups of the time-frequency plane using algebraic methods 
instead of complicated multi-index calculations. The understanding of these twisted 
group algebras is also intimately related to the symplectic structure of the time- 
frequency plane. Consequently we have devoted a part of this survey to describe 
harmonic analysis on the time-frequency plane. As our note is mostly addressed to 
applied mathematicians as well as engineers we did not go so far as to explain in 
which sense the Morita equivalence of the twisted group algebras of the subgroup 
and its adjoint subgroup. However, from a deeper view-point it can be seen as the 
very reason for the validity of duality theory of Gabor frames (see (30] for details). 
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